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In this paper we revise a perfect fluid FRW model with time-varying constants "but" taking 
into account the effects of a "c-variable" into the curvature tensor. We study the model under the 
following assumptions, div(T) = and div(T) ^ 0, and in each case the flat and the non-flat cases 
are studied. Once we have outlined the new field equations, it is showed in the flat case i.e. K = 0, 
that there is a non-trivial homothetic vector field i.e. that this case is self-similar. In this way, we 
find that there is only one symmetry, the scaling one, which induces the same solution that the 
obtained one in our previous model. At the same time we find that "constants" G and c must verify, 
as integration condition of the field equations, the relationship G/c 2 = const, in spite of that both 
"constants" vary. We also find that there is a narrow relationship between the equation of state and 
the behavior of the time functions G, c and the sign of A in such a way that these functions may be 
growing as well as decreasing functions on time t, while A may be a positive or negative decreasing 
function on time t. In the non-flat case it will be showed that there is not any symmetry. For the 
case div(T) 0, it will be studied again the flat and the non-flat cases. In order to solve this case 
" it is necessary to make some assumptions on the behavior of the time functions G, c and A. We also 

, find the flat case with div(T) — 0, is a particular solution of the general case div(T) ^ 0. 
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I. INTRODUCTION 

■ Since the pioneering work of Dirac ([!]), who proposed, motivated by the occurrence of large numbers in Universe, a 
theory with a time variable gravitational coupling constant G, cosmological models with variable G and nonvanishing 
and variable cosmological term have been intensively investigated in the physical literature (see for example [2]-[14]). 

Recently , the cosmological implications of a variable speed of light during the early evolution of the Universe have 
been considered. Varying speed of light (VSL) models proposed by Moffat ([15]) and Albrecht and Magueijo ([16]), in 
>) ■ which light was travelling faster in the early periods of the existence of the Universe, might solve the same problems as 
| inflation. Hence they could become a valuable alternative explanation of the dynamics and evolution of our Universe 

■ and provide an explanation for the problem of the variation of the physical "constants" . Einstein's field equations 
q ' (EFE) for Friedmann-Robertson-Walker (FRW) spacetime in the VSL theory have been solved by Barrow ([17]), 
q-i who also obtained the rate of variation of the speed of light required to solve the flatness and cosmological constant 

1 problems (see J. Magueijo ([18]) for a review of these theories). 
j^jQ, This model is formulated under the strong assumption that a c variable (where c stands for the speed of light) does 
r* ■ not introduce any corrections into the curvature tensor, furthermore, such formulation does not verify the covariance 
and the Lorentz invariance as well as the resulting field equations do not verify the Bianchi identities either (see 
Bassett et al [19]). 

Nevertheless, some authors (T. Harko and M. K. Mak [20] and P.P. Avclino and C.J. A. P. Martins [21]) have 
proposed a new generalization of General Relativity which also allows arbitrary changes in the speed of light, c, and 
the gravitational constant, G, but in such a way that variations in the speed of light introduce corrections to the 
curvature tensor in the Einstein equations in the cosmological frame. This new formulation is both covariant and 
Lorentz invariant and as we will show the resulting field equations (FE) verify the Bianchi identities. We will be able 
to obtain the energy conservation equation from the field equations as in the standard case. 

The purpose of this paper is to revise a FRW perfect fluid model with time varying constants (see [22] and [23]) but 
taking into account the effects of a c variable into the field equations. In our previous models ([22]- [23]) we worked 
under the assumption that a c— variable does not induce corrections into the curvature tensor and hence the classical 
Friedman equations remain valid. We will show that such effect is minimum but exists. In section 2, once we have 
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outlined the Einstein tensor, we check that it verifies that its covariant divergence vanishes. Consequently we impose 
the same condition to the right hand of the Einstein equation i.e. div (%2!G_jnij _|_ g^'A) = 0. In this way we set the new 
field equations as well as the Kretschmann scalars. As we will see, in order to integrate the resulting field equations 
we will need to consider some assumptions. In section 3, we will consider the condition div (T* J J = 0, and it will be 
studied two particular cases, the flat case, i.e. K = 0, and the non-flat case, K ^ 0. Since with the proposed method 
we are not able to solve the so called flat problem we need to study separately both cases (the flat and the non-flat 
cases) introducing such conditions as an assumption. 

In this approach, in the case K = 0, we will show that the model is self-similar since it is found a non-trivial 
homothetic vector field. In order to integrate the FE under the assumptions div(T) — and K = 0, we use the 
Lie group tactic which allows us to find a particular form of G and c for which our FE admit symmetries i.e. are 
intcgrable. As we will show under these assumptions the FE, only admit one symmetry, the scaling symmetry. This 
is the main difference with respect to our previous approach ([22]- [23]) where we did not consider the effects of a 
c-variable into the curvature tensor in such a way that the resulting FE admitted more symmetries. We also obtain 
as integration condition that the "constants" G and c must verify the relationship G/c 2 = const, in spite of that both 
"constants" vary. In this work we have found three solutions. The first one is very similar to the de-Sitter solution 
i.e. the energy density vanishes, the cosmological constant is a true constant while G and c follow an exponential law 
as the scale factor. The second and third solutions behave as the standard FRW where all the quantities follow power 
law with respect to time t, nevertheless the third solution is non-singular. In the second studied case, the non-flat case 
K 7^ 0, we will show that the FE do not admit any symmetry. Nevertheless we will try to find a particular solution 
imposing some restrictions. 

In section 4 we consider the possibility that div(T) ^ 0, and we will study again two particular cases, the flat case, 
i.e. K = 0, and the non-flat case, K / 0. The possibility that the covariant conservation condition div(T) = be 
relaxed has been advanced by Rastall ([24]), who pointed out that a non-zero divergence of the energy-momentum 
tensor has not been ruled out experimentally at all yet. In order to integrate the resulting FE we will need to make 
some assumptions (scaling assumptions) about the behavior of the "constants" G, c and A, which allow us to find 
particular solutions to the FE. This scaling assumptions work well in the flat (homothetic) case, as it is expected, 
but they seem very restrictive in the non-flat case. It is also showed that we can recover the solution obtained in the 
div(T) = as a particular solution of this model in such a way that G and c must verify the relationship G/c 2 = const. 
etc... as it is expected. 

In section 5 we end with a brief conclusions. 



II. THE MODEL 



We will use the field equations in the form: 
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where arbitrary variations in c and G will be allowed. We assume that variations in the speed of light introduce 
corrections to the curvature terms in the Einstein equations in the cosmological frame. In our model variations in the 
velocity of light are always allowed to contribute to the curvature terms. These contributions arc computed from the 
metric tensor in the usual way. The line element is defined by (we are following the O'Neill's notation [25]): 



ds 2 = -c(t) 2 dt 2 + f 2 (t) 



dr 2 
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and the energy momentum tensor is: 



Tij = {p + p) UiUj - pgij, 



(2) 



(3) 



where p and p satisfie the usual equation of state, p = mp in such a way that uj = const., usually to is taken such 
that, to € ( — 1, 1] , that is to say, our universe is modeled by a perfect fluid. The 4— velocity u % is defined as follows 
u l = (c- 1 ^), 0,0,0) such that u z u l = -1. 



The field equations are as follows: 
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where as we can see, the only difference with respect to the usual one is the factor —2 — 4- in eq. (4). Here a prime 
denotes differentiation with respect to time t. 

Since the divergence of [Rij ~ ^gijR) vanishes then we impose that the right hand of equation (1) has zero diver- 
gence too. Therefore applying the covariant divergence to the second member of the field equation we get: 

div (^-T ij + g ij A) = 0, (6) 



where we are considering that 22Jl j s a function on time t. Hence simplifying it yields: 



where H = f'/f. 

Therefore the new field equations are as follows: 

IE' -2-H + 3H 2 + ^ = -^fp + Ac 2 , (8) 
tt2 Kc 2 8ttG 1 A 2 
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As we have mentioned in the introduction these field equations are not new in the literature. They have been 
outlined by T. Harko and M. Mak ([20]) where they study a perfect fluid model and the influence of the time variation 
of the constants in the matter creation. Later they study Bianchi I and V model with time varying constants (taking 
into account the influence of a c-variable into the curvature tensor). Other author are P.P. Avelino and C.J. A. P. 
Martins ([21]) and H. Shojaic and M. Farhoudi, have obtained similar equations ([26]) since these authors consider G 
as a true constant. 

We would like to emphasize that deriving eq. (9) and substituting this result into eq. (8) it is obtained eq.(7) 
i.e. the covariant divergence of the right hand of our field equation in the same way as in the standard cosmological 
model. 

Curvature is described by the tensor field R^i- It is well known that if one uses the singular behavior of the tensor 
components or its derivates as a criterion for singularities, one gets into trouble since the singular behavior of the 
coordinates or the tetrad basis rather than the curvature tensor. In order to avoid this problem, one should examine 
the scalars formed out of the curvature. The invariants RiemS and RiccS (the Kretschmann scalars) are very useful 
for the study of the singular behavior, being these as follows: 



RiemS = R tl i rn R ljlrr \ (11) 



RiemS := —r 



12 

T 



c K 



(12) 



and 



RiccS := RijR 13 , 
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In order to try to solve eqs. (8)-(10) we need to make some hypotheses on the behavior of the quantities. 



(13) 



III. SOLUTION I. divT = 0. 



In this first solution we make the following assumption 

divT = (14) 

and we will consider two subclasses, the flat and the non-flat cases. 

Our tactic consists in studying the field equations through the Lie method (see [27] — [28] — [29]). In particular we 
seek the forms of G and c for which our field equations admit symmetries i.e. are integrable. 

We will consider the following assumption divT = 0, which transforms eq.(10) into these two new equations: 

p' + Z{oj + l)pH = 0, (15) 
G' „c' A'c 4 

G- A ---^G P (16) 

In order to use the Lie method, we rewrite the field equations as follows. From (8 and 9) we obtain 
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From equation (15), we can obtain 
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where A = 12tt (u> + 1) > 0, Vw. Following this tactic we try to make the smallest hypothesis number and to obtain 
the exact behavior of the "constants" G, c and A. Following other tactics we are obliged to make assumptions that 
could be unphysical. 



A. Flat case. Self-similar approach. 

As have been pointed out by Carr and Coley ([30]), the existence of self- similar solutions (Barenblatt and Zeldovich 
([31])) is related to conservation laws and to the invariance of the problem with respect to the group of similarity 
transformations of quantities with independent dimensions. This can be characterized within general relativity by 
the existence of a homothctic vector field and for this reason one must distinguish between geometrical and physical 
self-similarity Geometrical similarity is a property of the spacctime metric, whereas physical similarity is a property 
of the matter fields (our case). In the case of perfect fluid solutions admitting a homothctic vector, geometrical 
self-similarity implies physical self-similarity. 

As we show in this section as well as in previous works, the assumption of self-similarity reduces the mathematical 
complexity of the governing differential equations. This makes such solutions easier to study mathematically. Indeed 
self-similarity in the broadest Lie sense refers to an invariance which allows such a reduction. 

Perfect fluid space-times admitting a homothctic vector within general relativity have been studied by Eardlcy 
([ ]). In such space-times, all physical transformations occur according to their respective dimensions, in such a 
way that geometric and physical self-similarity coincide. It is said that these space-times admit a transitive similarity 
group and space-times admitting a non-trivial similarity group are called self-similar. Our model i.e. a flat FRW 
model with a perfect fluid stress-energy tensor has this property and as already have been pointed out by Wainwright 
([33]), this model has a power law solution. 



Under the action of a similarity group, each physical quantity <fi transforms according to its dimension q under 
the scale transformation. For space-times with a transitive similarity group, dimcnsionless quantities are therefore 
spacetime constants. This implies that the ratio of the pressure of the energy density is constant so that the only 
possible equation of state is the usual one in cosmology i.e. p = up, where u is a constant. In the same way, the 
existence of homothetic vector implies the existence of conserved quantities. 

In the first place we would like to emphasize that under the hypothesis K = 0, we have found that the space-time 
(M,g) is self-similar since we have been able to calculate a homothetic vector field, Xh £ X(M), that is to say 

L x g = 2g, (21) 

where 

Xh _ (if^tft) a + ( x _ (Mg±^) n) (* + yO, + *) . (22) 

This kind of space-times have been studied by (Eardly ([32]), Wainwright ([33]), K. Rosquits and R. Jantzen ([34]) 
and for reviews see Carr and Coley ([30]) and Duggal et al ([35])). 

As we have mentioned above we are only interested in the case K = for this reason eq. (20) yields 



p"= P - + /-+A%p\ (23) 



Now, we apply the standard Lie procedure to this equation. A vector field X = £(i, p)d t + rj(t, p)d p , is a symmetry 
of (23) iff 

Hp + pt PP = 0, (24) 
rjc - 2c'p 2 £ p - cpri p + cp 2 ij pp - 2cp 2 £ tp = 0, (25) 
-cc'p^t - 3p 3 AG£ p + (c' 2 p - cc"p) £ + 2c 2 p Vtp - c 2 p^ tt - 2c 2 i lt = 0, (26) 
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Solving (24-27), we find that 

£ = at + b, 77 = -lap, (28) 



subject to the following constrains, from eq. (26): 



and from eq. (27) 
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where a and b arc numerical constants. 

We would like to emphasize that this is the main difference with respect to our previous work ([22]) where we found 
that the model admitted more symmetries. Note that in this paper we are considering that a c-variable introduces 
corrections into the curvature tensor, this possibility brings us to obtain a new eq. (23) which contains an extra term 

{f'l^j w ith respect to the employed one in our previous paper ([22]). 

In this situation we have found that this model only admits two symmetries X\ = d t i.e. a movement and 
X2 = tdt — 2pd p , i.e. a scaling symmetry as it is expected in this kind of models (self-similar model). In such a way 
that [Xtj-X^] = X\ i.e. they form a L2 Lie-algebra, see ([28]). 

From (30) we can see that 
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(where we will assume that B > 0) that is to say, that both constants vary but in such a way that the relationship ^ 
remains constant for any t and independently of the any value of the constants a and b. For this reason our equation 
(23) may be rewritten as follows 

p" = P —+p'- + A P \ (32) 
p c 

where A = 12ir (uj + l) 2 B = const > 0, Vw, (u> ^ —1), in such a way that all the restrictions come from eq. (29). 

The knowledge of one symmetry X might suggest the form of a particular solution as an invariant of the operator 
X i.e. as solution of 

% = "A (33) 

this particular solution is known as an invariant solution (generalization of similarity solution) furthermore an invariant 
solution is in fact a particular singular solution. 

In order to solve (32), we consider the following cases. 

1. Case I. 

Taking a = 0, b ^ 0, we get 



c" = — , c(t)=K 2 e K i\ (34) 

c 

where the {Ki) 2 i=l are integration constants. In this way we find that from eq. (30) G behaves as: 

G = /i|e 2Klt . (35) 

Therefore the equation (32) yields: 

p" = ^+p'K l +Ap 2 , (36) 
P 

where X\ = dt- The use of the canonical variables brings us to an Abel ode: 

?/ = -Ax 2 y 3 -K iy 2 -l, (37) 
x 

where x = p and y = 1/p', since we cannot solve it, then we try to obtain a solution through the invariants 
This method brings us to the following relationship 

dt dp . . 

— = — ==>• p ~ const := poj (38) 

which seems not to be physical. This particular solution must satisfy eq. (32) which means that p — 0. With this 
solution we go back to eq. (9) 

3iJ 2 = Ac 2 , (39) 
and from eq. (16) it is obtained the behavior of "constant" A 
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with 



RiemS := 2AK$, RiccS := 36Kf. (42) 

This symmetry has brought us to obtain a super de Sitter solution, with a energy density vanishing and with a 
scale factor growing like a power of exponential functions, while G and c follow an exponential functions that depends 
of the constant K\. If for example we fix K\ = 1, then there is a sudden singularity in a finite time, in this case 
both "constants" G and c are growing functions on time t. While if K\ = — 1, the scale factor reaches an asymptotic 
behavior with G and c decreasing on time t. In any case A = const. > 0, i.e. is a true positive constant. 



2. Case II. 



Taking b = 0, a ^ 0, we get that the infinitesimal X is X\ = td t — 2pd p , which is precisely the generator of the 
scaling symmetries. Therefore the invariant solution will be the same than the obtained one with the dimensional 
method. 

With these values of a and b equation (29) yields 

c" = — -y, => c(t) = K 2 t K \ (43) 
c t 

as we expected, c(t) follows a power-law solution (self-similar solution), where K\ and K 2 are numerical constants, 
Ki, K 2 0. We would like to emphasize that with this behavior for c the homothetic vector field behaves as: 

XH = l^ dt + { l -l^Tl) {xdx + yd y + zdz) (44) 

as in the standard model (see for example Eardley ([32]) and Wainwrith ([33])). 
In this way G behaves as: 



G = K 2 t 1 . (45) 



Therefore equation (32) yields: 



p"= P -+p'^+Ap\ (46) 
P t 



we expect that this ode admits a scaling solution (i.e. that the energy density follows a power law solution) as it is 
expected in this kind of models. 

The canonical variables bring us to obtain an Abel ode 

y' = x (2 + 2K X - Ax) y 3 - (1 + I<i)y 2 - l , (47) 



x 



where x = pt 2 and y = t 2^ p ,\ +2 p^ , and which solution is unknown 
The solution obtained through invariants is: 

dt dp o 2(l + ATi) 

finding that K\ > — 1 from physical considerations i.e. p > iff Ki > — 1. This is the kind of solution expected in a 
self-similar model (see for example Wainwrith ([33]) and Jantzen ([34])). 
From eq. (16) it is obtained the behavior of "constant" A 

c' WirBKx , . 
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as it is observed if K\ > then A is a negative decreasing function on time t, while if K\ < then A is a decreasing 
function on time t, but with the restriction K\ 6 ( — 1,0). We would like to emphasize that the self-similar relationship 



is trivially verified as it is expected in this kind of solutions. 

Now, we will calculate the scale factor /. In order to do that we may follow two ways, in the first one, from eq. (15) 
we find that 



/ = /ot 2/3(.+l) ) (51) 
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where /o = ^ 2(1 + K ) ) ' °^ servm S a S am that necessarily K\ > — 1. 

If we follow our second way it is found that from eq. (9 with K = 0) 
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^/Hq = —. — — , as it is expected from eq. (51). Therefore this new way does not add any more information. 
In this way we found that 

H=y/H t-\ and q = | (l) - 1 = - 1. (53) 

therefore q < iff u) £ (—1, — |) , note that if u) < — 1 (phantom cosmology), then f(t) is a decreasing function on 
time t. 

If we make the assumption (scaling symmetry) on the scale factor 

fact (54) 
in such a way that equating this expression with eq. (51) then 

< 55 > 

finding that K 1 > uo £ ( — 1, — |) . Hence, if K\ > w £ (— 1, —5), G and c are growing functions on 
time t, while A < 0, i.e. is a negative decreasing function on time t. In other case, if K\ < 0, G and c are decreasing 
functions on time t, while A > 0, i.e. is a decreasing function on time t. 
The Kretschmann scalars behave as: 

RiemS ~ K ltil+Ki) ' RlccS ~ K^i+W) ' (56) 

finding in this way that if K\ < — 1 (forbidden possibility) then both scalars tend to zero while if K\ > —1, then both 
scalars tend to infinity i.e. there is a true singularity. 

In the first place we would like to emphasize that this is the solution that we have obtained in our previous works 
using D.A. (see [23] and [36]). This is due to two reasons. The first one, because D.A. is a powerful tool that may 
be used even when the FE are not well outlined. The second one, because in the previous works we were using only 
three of the FE, ignoring eq. (8) and using eqs. (9,15 and 16). Therefore the solution found in those papers work 
well and we refer to them for all the physical considerations. 

In this approach we are supposing that q < i.e. the universe accelerates due to a equation of state ui £ (—1, —3) , 
but as it has been pointed out by R.G.Vishwakarma (see [37]) the acceleration of the universe may be explained 
through different mechanisms in such a way that q < while ui £ [0, 1] . In this way we would like to stress that 
to solve some of the cosmological problems that present the standard model, we have found that K\ > —1, with 

A'i = — - ^ in such a way that if K\ >0iffw£(— 1, — 5) then q < and G and c are growing functions on time 

i, while A is a negative decreasing function on time t. With K\ < 0, q > and G and c are decreasing functions on 
time t, while A is a positive decreasing function on time t, i.e. there is a narrow relationship between the behavior 
of G and c and the sign of A controlled by oj. 



3. Case III. 



This case is a generalization of the above case, simply, we will avoid the singular case. Taking a, b ^ we get 

„'2 n J K, 
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obtaining in this case that G behaves as: 



c(t) = K 2 (at + b)~ , (57) 



hence equation (23) yields: 



G = Ki(at + b)~, (58) 
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The canonical variables bring us to an Abel ode 

y' = x (2a 2 + 2aK x - Ax) y 3 - (a + K x )y 2 - - (60) 



where x — p(at + b) and y = ^ at+b ^ ^ p ,^ at+b ^ +2a p) i that is to say, we obtain a very complicate ode which at this 
time we do not know how to solve. Nevertheless, as we have pointed out in ([23]) the general solution of this kind of 
equations are unphysical, for this reason it is sufficient consider particular solutions obtained through the invariants. 
The solution obtained through invariants is: 



dt _ dp _ 2a(a + K x ) 

T ~ ~V^ P ~ A(at + bf 
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where it is observed that K\ > —a (but if you choose a = b = 1, then we have the same result as before i.e. K\ > —1). 
Therefore the scale factor behaves as: 

/ 2 s 1/3(oj+1) 

[ A U A (at + b) \ ,x2/3(w+i) / RO x 

/= (^ + J ' / = /o(a* + 6) ( 62 ) 

which is very similar to the last result (see eq. (51) but in this occasion this solution is non-singular). 
We end finding the behavior of A, as in previous cases taking into account eq. (16) it is obtained: 

. , „ c' 167r BaK\ 
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as we can see in this case if K\ > then A is a negative decreasing function on time t. 
The Krctschmann scalars behaves as: 

RimS « .,- , „ , , RiccS « ... , „ , , (64) 

X|(at + &) 4(1+ifl) tf| (o* + 6) 4(1+Xl) 

showing a non-singular state when t runs to zero. 

This solution is very similar to the previous one except that this solution is non-singular. In fact, the above solution 
is a particular solution of this one. 



B. The non-flat case. 



In this subsection we go next to study the particular case K ^ 0. One of the drawbacks of the above approach is 
that we need to make the assumption K = i.e. our approach is unable of solving the so-called flatness problem. 
In order to research if it is possible to solve such problem we go next to study eq. (20) through the Lie method, 
seeking symmetries that allow us to obtain any solution in closed from. But as we will see eq. (20) does not admit 
any symmetry in such a way that in order to obtain a particular solution we will impose a concrete symmetry, but 



this is precisely the method that we are trying to avoid, to make assumptions or at least to make the minor number 
of assumptions or to make assumptions under any physical or mathematical (symmetries) well founded reasons. We 
only explore one case. 

Therefore the equation under study is: 
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P" = — +p'--2—p + A-p\ (65) 



but as 



and hence 



tivT = 0^ p = A u f^ +1 \=> f= [-£-)""'" ,=> f = (p)*&m , (66) 



p — +p ' C --2c 2 p a + A% P \ (67) 
p c 



where a = g^jqrj-y, and for simplicity we have adopted the case K = 1. 
The Lie group method brings us to obtain the following system of pdes 



p£pp + = 0, (68) 
-p-\ - 2-£ p + r} pp - 2£ tp + p- 2 v = 0, (69) 



r> \ 2 r" \ f C 

C \ C \ ^ „ 9 / „ n 2 , (_r 



- U + 3p 2 2c z p a ~ 2 - A— £ p - 2 P - L r lt - + 2 Vtp £ t =0, (70) 



C ' 2 ( a G o„2„a-2\„ , r,2(r,2a-2 ' @ 



Vtt - -Vt+P z [A- ~ 2cV^ I Vp + 2p< [2c<p a -< - A— ) £ 



which has no solution, that is to say, eq. (67) does not admit any symmetry, for this reason we will need to follow 
other approaches. 

For example, if we "impose" any particular symmetry X, maybe we may found some restrictions for the behavior 
of the quantities G, c and p. We will explore such possibility. 



1. Case I. 

In this case, we choose (£ = l,r] = 0) ,i.e. X = d t , in such a way that from eq.(70) it is obtained the following 
restriction 

r' 2 

c" = — => c(t)=K 2 e K >\ (72) 
c 

and from eq.(71) it is obtained the following one 



hence 



^ A^ + 1 ) Kl = %> (?4) 



as we can see from eq. (73), we cannot obtain the condition G = Be 2 (as in the flat solution) since such condition 
means that 

„ o c 4 c G c 
Ap AG-=G- 2 - = °^ P = °> (75) 

that is to say, the energy density vanishes. 

In order to find a particular solution to eq. (74) we impose the condition a = 2 (as mathematical condition) which 
means that ui = — f <C — 1, although such possibility is very restrictive (and maybe unphysical, the ultra phantom 
equation of state). In this way it is found that 

G' = IKxG+^^f- ,=>G(t) = C ie 2tK i + 4^ 2 4 e 4(t/fl) , (76) 

where C± is an integration constant, therefore eq. (67) yields 

p" = P ^+ P 'K 1+ [A^p\ (77) 

since this equation has not an explicit (analytical) solution, note that we have obtained the same eq. as in section 3.1 
case I eq. (36), we find again that a particular solution is p = 0, therefore from the field equation (9) 

H 2 + = y*^ 1 *. (78) 

and hence in this case we have found that a solution is: 

f = —ff= o + 7 CX P ~, t? (cxp(ifit) + Ci) cxp — (cxp(Ait) - Ci) , (79) 



/So V 3 6 1 \ 3 Ki K CK ' ' ) r \ 3 K 



( G' c'\ A'c 4 
= — - 4- )p = -—— => A = A = const. > 0. (80) 
\ G c J 87rG 

This solution looks very unphysical or at least very restrictive, a vanishing energy density and a constant cosmolog- 
ical constant and is very similar to the obtained one in section 3.1. case I the super de-Sitter solution, with the same 
behavior and the same restriction for the numerical constant K\, i.e. sudden singularities or asymptotic behavior 
depending of the sign of K\. In this case, the employed method does not allow us to obtain the behavior of G, c and 
A since the equation under study does not admit any symmetry. 



IV. SOLUTION II. divT / 0. 

In the previous section we have made the assumption that div(T) = i.e. the divergence of the energy-momentum 
tensor vanishes. Nevertheless we have obtained as general conservation equation eq. (10) i.e. 

/ G' c'\ A'c 4 , . 

P' + 3(. + l)^ = -(--4-J^_ , (81) 

i.e. div(T) = p' + 3 (u) + 1) pH and we have assumed a particular case (with perfect mathematical sense) 
div(T) = 0. In this section we will study the general case div (T) ^ 0. 

The possibility that cosmological and physical considerations may require that the covariant conservation condition 
div(T) = be relaxed has been advanced by Rastall ([24]), who pointed out that a non-zero divergence of the energy- 
momentum tensor has as yet not been ruled out experimentally at all. In Rastall's theory ([24]), the divergence of T is 
assumed to be proportional to the gradient of the scalar curvature S, div(T) = Xgrad(S) , where a A = constant, and 
in fact the modified field equations are equivalent to standard general relativity with an additional variable A term. 
We refer to the reader to the Harko and Mark work ([20]) to see a matter creation and thermodynamical approach in 
this context. 



A. The flat case. 



In this case we are going to consider the field equations "but" without the condition div(T) = 0, i.e. the field 
equations will be 

2H' -2-H + 3H 2 = -^-p + Ac 2 , (82) 
c c l 



3H 2 = -^ P + Ac 2 , (83) 
c l 

✓ + S fr, + l),tf_-(§-4),-|g, (84) 

but as we can see we have 2 equations with 5 unknowns, therefore it is necessary to make some assumptions, for 
example in the behavior of G, c and A as follows: 

G = G a H a , c = c H b , A = A c- 2 H 2 , (85) 

where Gq,cq and Ao are dimensional constants and a, 6 £ R, without any restriction i.e. we do not need to assume 
any concrete sing or value for these numerical constants. Furthermore, we must stress that the conditions K = 0, 
together to power law assumptions bring us to a scaling solution. 

Taking into account these assumptions and form eq. (83) we obtain p 



P 



3-A 



ff 2 < 1+b )- a , i.e. p = Po H^ 1+b '>- a , (86) 



d 

where do = 87tGo/cq, and taking this relationship into eq. (84), we obtain the following ode in quadrature 

a^ + 3(u; + l)H = -a^+4b^-K^, (87) 

and therefore 

[a + a - 46 + Kj — = -3 (u + 1) , 

where a = 2(1 + 6) — a, K = (2(1 — 6)coA /87rGo/Oo) , and therefore 
H' _ (w + l)(3-Ao) 



H 2 2(1-6) 



H = h nt-\ f = foF° n , (89) 



where n = - — ^.j - — — — , and ho = const > 0, is an integration constant and we impose that 6 7^ l,Ao 7^ 3 and 

w ^ -1. As we can see, there are some restrictions, for example, if lj G (—1,1] then 6, Ao must verify at the same time 
that 6 < 1 and Ao < 3 or 6 > 1 and Ao > 3. Now if to < —1 then it should exist a combination between the signs of 
6, Ao such that ?? > 0, in other way, the radius of the Universe decreases. It is observed that the behavior of the scale 
factor does not depend of the constant a, i.e. of the behavior of the gravitational "constant" . 

Once we have obtained the behavior of / i.e. of H we go next to complete our calculations of the rest of the 
quantities i.e. 



with the restrictions 



poH 2(i +b) -a ^ p = Po (h nt- 1 ) 2(1+b) - a »t a - 2 ^+ b \ (90) 



po (h nf {1+b) a > 0, and a - 2(1 + 6) < 0, (91) 



i.e. we are assuming that the energy density is a positive decreasing function on time t. 

It is observed too that if 26 = a, then we obtain the particular solution p « t~ 2 , as well as the realtionship 
G/c 2 = const, as in the above cases, i.e. the obtained solution under the assumption div(T) = 0. In this way we can 
se that the div(T) = case is a particular solution of the div(T) ^ case, as one may expected, but we have not any 
physical or mathematical (symmetry or intcgrability condition) reason to assume such relationship. 



We end calculating the behavior of the "constants" G, c and A i.e. 



G = GoKn a t- a , c = c {h n) b t~ b 1 A = A c^ 2 (/i n) 2(b_1 ' 1 t 2 ' fc_1) , (92) 

in such a way that A will be a decreasing function on time iff b < 1. 
In this way we found that 

i d ( 1 \ 1 2(1-6) 

H = h nt-\ and q=-(-)-l = - 1= v y - 1, 93 

at \H J h§n h (uj + 1) (3 — A ) 

therefore g < iff |2 (1 - 6)| < |/i (w + 1) (3 - A )| . 
The Kretschmann scalars behave as: 

RiemS « . . — rr, RiccS ~ , ... — rr, (94) 

c 4 f 4(l-6)' C^ 1 -^ ' 

finding in this way that there is a true singularity. 

In this case we have found a scaling solution as it is expected for this kind of models (K = 0). As we can see the 
scale factor / docs not depend of G only of c and A. The energy density p depends of G and c. Nevertheless we have 
not been able to find better restrictions for the introduced (ad hoc) numerical constants a and 6 such that they give 
us more information about the behavior of the time functions G, c and A 



B. The non-flat case. 



As in the above section we will study the case K ^ separately. For this purpose the FE are now: 

2H' -2-H + 3H 2 + ^ = -^fp + Ac 2 , (95) 

+■ ^Ac 2 , (96) 

c'\ A'c 4 , . 

" 4 - P- 5-77. ( 97 ) 
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but as we can see we have 2 equations with 5 unknowns, therefore it is necessary to make some assumptions, for 
example in the behavior of G, c and A as follows: 

G = G a H a , c = c Q H b , A = A QC - 2 H 2 , (98) 

where Go, cq and Ao are numerical constants. We must stress that in this occasion these hypotheses could be unphysical 
since K ^ and we are imposing a scaling behavior typical of the flat case. 
Taking into account these assumptions and form eq. (96) we obtain p 

p = (^) + U . p = ^««>- + ^, (99) 

where do — 87tGo/cq, and taking this relationship into eq. (97), we obtain the following second order for / 

(f') 2b (f) 2 

f' = - D i)Um +D2l r> (100) 

where 

_ (3a; + 1) Kc 2 f (tt + 1) (3-A ) \ 

Dl ~ 2 (1-6)' ^-{^ 2 (1-6) J' (101) 

In the first place we may note that (6 ^ 1) and that if u = —1/3 D\ = 0, independently of the value of constant 
K, in this case we are mainly interested in the K ^ case. If K = 0, or D\ = then we obtain again the solution 
already obtained in the latter (last) case. 



Calculation of Eq. (100). making the following change of variables it is obtained the first order ode 

(x = f, y = j) =► y' = D 1 x-^y*-* b -D 2 x- l y, (102) 



and which solution is: 



therefore 



1 (26£> 2 - 2L> 2 - 26) ' 



f , = ( r f 2D 2 (i- b ) , 2(6 -!)£>! 2b \ * 

^ V (2bD 2 -2D 2 -2b) J 



(104) 



and hence 



'=/ ( c -- w+ (CTi) -i" A+ ' i ' <io5) 

Since we have not obtain information about the behavior of / then we try to find a particular solution. For this 
purpose, we observe that eq. (100) admits the following symmetries: 

X 1 = 8 t , X 2 = td t + (l-b)fd f , (106) 

where we would emphasize that X 2 is a scaling symmetry, maybe induced by the hypotheses about the behavior of 
G,c and A, (scaling relationships). 

The invariant solution (particular solution) that induces X 2 is the following one: 

dt df f (1 _ 6) 



t (1-6)/ 

with 



f = fot (1 - b) , (107) 



1M6- 1) l D 2 {b-l)-b _ A-cg(l + 3a;) / (3a; + 1 - A (q; + 1)) _ 
/0 D 2 (6-1)-6Y £>i(6-l) ~ (3cj + l-A (w + l))y A'cg (1 + 3lu) C ° nS ^ {WS) 

and we necessarily impose that b < 1 and 6^0, this means that if 6 — > 1 them / — > /o. As we can observe K ^ 0, 
and we must be careful with the signs since if AT = —1 then it should be satisfied the relationship 3ui + 1 < Aq(uj + 1). 
In this way we found that 

H=(l- b) t-\ and ,= 1(^-1 = ^-1 = ^, (109) 

therefore g < iff 6 < 0, which implies that c is a growing function on time t. 
If we make the assumption (scaling-symmetry) on the scale factor 

/wet /«^-«, (110) 

The Kretschmann scalars behave as: 

ffiemff « ^^i) ' fliocgw g^t-i) ' ( U1 ) 

finding in this way that there is a true singularity. 

Calculation of p, and the possible restrictions for the constants a and 6. 

p = po (1 - 6) 2h - Q+2 t-a(W-i) + (i _ 4 -2 ; (n2) 

/o 

and 

G«G *~ a , c«c i" & , A«A^ W) . (113) 

As we have seen, we have only been able to obtain a particular solution. Imposing so many hypotheses, we lose 
information and we do not know how to recover it in such a way that we are not able to know better the behavior of 
the time functions G, c and A. 



V. CONCLUSIONS. 



In this paper we have studied a perfect fluid FRW model with time-varying constants "but" taking into account 
the possible effects of a c— variable into the curvature tensor. In this way, as other authors have already pointed out, 
such effects are minimum in the field equations but they exists and are very restrictive. Under the made hypotheses, 
we have seen that the Einstein tensor has covariant divergence zero, in this way we have imposed that the right hand 
of the field equations has a vanishing divergence too i.e. 



In this way we have obtained the set of the new FE, emphasizing the fact that we can recover eq. (114) from the 
field equations as in the standard case i.e. deriving one of them and simplifying with the other one. 

In other to solve the resulting FE we have considered the following cases. In the first case we have imposed the 
condition div(T) = 0, as a particular case of eq. (114) and we have studied the flat and non-flat cases. We have 
needed to make such distinction because with the employed method we are not able of solving the so called flatness 
problem, for this reason we have needed to study separately then. 

The flat case under this hypothesis is the already studied one in our previous works (see [22]) and under this 
new considerations i.e. taking into account the effects of a c-variable into the curvature tensor, we have shown that 
the scaling solutions the only one while in ([22]), without this new assumption, we obtained more solutions, i.e. we 
obtained other solutions apart from the scaling solution. We would like to emphasize that it has been obtained as 
integration condition that the "constants" G and c must verify the relationship G/c 2 — cosnt. in spite of the fact that 
both "constants" vary. This result is in agreement with our scaling solution obtained in our previous work (see [22]) 
but in those works we needed to make such relationship as assumption. 

We have also shown that this model is self-similar since we have been able of obtaining a non-trivial homothetic 
vector field. This result is agreement with the obtained scaling solution as it is well known. 

With the obtained solution, if we want to solve the acceleration of the universe i.e. q < 0, then the time function 
G and c are growing function on time t, while A is a negative decreasing function on time t, in this case the equation 
of state belong to the interval: uj £ (—1, —1/3). In other cases G and c are decreasing function while A is a positive 
decreasing function and in this case the equation of state belongs to the interval: u £ [—1/3, 1] . 

The non-flat case does not admit any symmetry and the particular studied case looks very unphysical. 

The second class of studied models verifies the condition div(T) ^ 0, i.e. without imposing any restriction to eq. 
(114) and we have studied again the flat and non-flat cases. To solve these cases we have needed to make scaling 
assumptions on the behavior of the time functions G, c and A. These assumptions work well in the flat case (self- 
similar case) but in the non-flat case seem very restrictive. In the flat case, it is found for the made assumptions, 
that there is a relationship between the numerical constants that determine the behavior of G and c that brings us to 
obtain again the particular case div(T) = 0, i.e. such case could be seen as a particular solution of the more general 
case div(T) ^ 0. 



We wish to acknowledge to Javier Aceves his translation into English of this paper. 
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